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Abstract 

In this paper, the operation of a Licensed Shared Access (LSA) system is investigated, considering 
downlink communication. The system comprises of a Multiple-Input-Single-Output (MISO) incumbent 
transmitter (TX) - receiver (RX) pair, which offers a spectrum sharing opportunity to a MISO licensee 
TX-RX pair. Our main contribution is the design of a coordinated transmission scheme, inspired by the 
underlay Cognitive Radio (CR) approach, with the aim of maximizing the average rate of the licensee, 
subject to an average rate constraint for the incumbent. In contrast to most prior works on underlay 
CR, the coordination of the two TXs takes place under a realistic Channel State Information (CSI) 
scenario, where each TX has sole access to the instantaneous direct channel of its served terminal. Such 
a CSI knowledge setting brings about a formulation based on the theory of Team Decisions, whereby 
the TXs aim at optimizing a common objective given the same constraint set, on the basis of individual 
channel information. Consequently, a novel set of applicable precoding schemes is proposed. Relying on 
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statistical coordination criteria, the two TXs cooperate in the lack of any instantaneous CSI exchange. We 
verify by simulations that our novel coordinated precoding scheme outperforms the standard underlay 
CR approach. 


Index Terms 

Spectrum sharing, coordination, precoding, local CSI, QoS 

I. Introduction 

The utilization of the radio speetrum is internationally regulated by governments, with the 
aim of providing wireless eommunieation serviees that ean be effieiently proteeted from harmful 
interferenee. Nevertheless, the tremendous spread of wireless serviees has given rise to a great 
need for bandwidth, whieh eannot be satisfied by an exelusivity of speetral alloeation. On the 
other hand, spectrum scarcity, i.e., the phenomenon where the radio speetrum is beeoming 
vastly underutilized, has been diseussed in reports, sueh as the one published by the Federal 
Communieations Commission (FCC) in 2002 [[I]|. The spotted under-utilization has given rise, 
in its turn, to the notion of Cognitive Radio (CR), whieh has been suggested as a promising 
teehnology in view of inereasing wireless speetral effieieney by exploiting the existing spectrum 
holes in time, frequeney or spaee Il2l. If3l. 

Foeusing on the underlay CR approaeh, a primary network allows the simultaneous use of 
its speetral resourees by a new-eoming (unlieensed) seeondary network, given the eondition 
that the latter will utilize the available resourees in a way that the interferenee ereated by a 
seeondary transmitter (TX) towards a primary reeeiver (RX) is below a threshold predefined 
by the primary network [@1, [|5||. Under sueh a setup, effieient sehemes have been proposed, 
emphasizing on multiple-antenna settings, with the aim of maximizing the information rate 
of the seeondary system, subjeet to given eonstraints over the harmful interferenee suffered 
by primary terminals [|^- lfT0l . Importantly, in praetiee, the ability of the seeondary TX to 
aequire global, multi-user Channel State Information (CSI) is unrealistie. At this point, one 
would suggest an exehange protoeol of ehannel information between TXs. However, sueh a 
protoeol would introduee exehange imperfeetions of ehannel estimates as well as delays, henee, 
leading to degraded performanee of the designed preeoding solution. Consequently, an extensive 
literature has emphasized on designing transmission sehemes, that are robust to imperfeet CSI or 
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merely requiring loeal ehannel knowledge [See [frill and referenees therein]. In addition, iterative 
sehemes, based on game theory, have been also investigated as a means of avoiding the need 
for global multi-user CSI exehange, with respeet to speetrum sharing seenarios l[T^ - l[T4ll . 

Yet, standard speetrum sharing approaehes for underlay CR systems, merely foeus on designing 
a transmission seheme for the seeondary TX, sueh that the interferenee reeeived by a primary 
terminal would not overeome a eertain threshold, beyond whieh a seeure primary eonneetion 
eannot be established. However, standardization bodies have lately foeused on the design of 
Authorized or Lieensed Shared Aeeess systems (termed as ASA and LSA) lITSll . ifT^ . The key 
differenee between the latter systems and underlay CR systems is that the ineumbents (equivalent 
to primary nodes in a CR system) ean share the speetrum with the lieensees (the lieensed 
equivalent of seeondary nodes in a CR system), provided that Quality-of-Serviee (QoS) metries, 
that have been negotiated prior to lieensing, are satisfied for all involved entities. Motivated by 
this new framework, it is evident that a major drawbaek of standard, interferenee temperature- 
based underlay CR systems eonsists in the laek of eoordination between the primary and the 
seeondary systems. As a result, in sueh a way, the primary system tends to overspend its available 
resourees, leading to poor throughput performanee at the seeondary side. 

Given this situation, in this work, we propose the design of a eoordination seheme for the 
two TXs, based on eommonly available, slow-varying statistieal information. The goal of sueh a 
design will be the maximization of the (average) information rate of the lieensee system, given 
that the aehievable (average) rate of the ineumbent system lies above a eertain threshold. More 
preeisely, aeeording to this seheme, the ineumbent TX exploits its loeally available CSI and 
the statistieal (eovarianee) information of the global multi-user ehannel to eoordinate with the 
lieensee TX, in order to guarantee the desired average throughput for its assigned terminal. Eaeh 
TX has aeeess to the instantaneous direet links of its assigned users (as well as their statisties), 
whereas, the interferenee eross-links are just statistieally known at both TXs. The availability of 
different estimates of the global downlink ehannel falls within the paradigm of Team Decision 
theory, beeause both transmitters (ineumbent and lieensee) are aetively engaged in eooperation, 
while being eonstrained by the loeality of the available instantaneous CSI llT7l - ll2TI . 

In [|22l . a similar seenario was investigated for a single-user underlay CR setup, in the existenee 
of spatially uneorrelated direet ehannel links, while, in the present work, we foeus on the perfor¬ 
manee of an extended set of applieable joint preeoding sehemes, with the assumption of eorrelated 
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Rayleigh fading for all the involved MISO channels. More particularly, our contributions are the 
following: 

• We design a low complexity, statistically coordinated precoding scheme for a MISO spec¬ 
trum sharing system, which can be applicable to a shared spectrum access system (ASA 
or LSA). The goal of this design is the maximization of the throughput of the licensee 
network, in terms of the achievable average rate, given a QoS constraint on the average 
rate of the incumbent. 

• In order to design this joint precoding scheme, we derive lower bounds for the average rate 
of each RX, in closed form, both when a precoding solution based on Matched Filtering 
(MF) and statistical Zero-Forcing (sZF), is applied. This way, since the original optimization 
problem is hard to solve, because of the existence of combined instantaneous and statistical 
CSI at each of the two TXs, we focus on an approximated version of it. 

• Focusing on the approximated version of the optimization problem, we design a set of 
applicable, joint precoding solutions, the elements of which are joint transmission schemes 
based on MF or sZF-based beamforming solutions. 

• For each of the joint transmission schemes that belong to the referred set, a power policy 
coordination criterion, which is based on the statistics of the global downlink channel, is 
applied, with the aim of finding the average transmit power levels at each TX. These power 
levels are such that the average rate criterion on the incumbent RX is satisfied. Then, the 
joint beamforming and power allocation scheme that maximizes the average throughput of 
the licensee RX is selected for transmission. Since such a decision relies on commonly 
available statistical information, it is taken by both TXs coherently. 

• The novel, coordinated precoding scheme is numerically evaluated in comparison to the 
standard interference temperature-based underlay CR approach. It is shown that our scheme 
outperforms the standard underlay CR one in a range of system scenarios, which makes 
our designed precoding policy appealing for LSA systems, since higher average throughput 
is achieved for the (candidate) licensee system. 

Throughout the paper, the following notations are adopted: all boldface letters indicate vectors 
(lower case) or matrices (upper case). A^, tr(A) and [A]^^ denote the Hermitian transpose 
of matrix A, its trace, and its (m,n)-th entry, respectively, whereas Aj(A) stands for its j- 
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th eigenvalue. Also, diag (ai,..., «„) symbolizes a diagonal matrix, the elements of whieh are 
ai,... ,an- Additionally, E[-] symbolizes the expeetation operator and || ■ || denotes the Euelidean 
norm, while 0„ denotes the all-zero veetor of dimension n. The identity matrix of dimension 
n X n is denoted by I„, while i denotes the eomplementary index of i, when the eardinality of 
the eonsidered set is equal to two, i.e., i = i mod 2 + 1. The phase angle between two veetors 
a and b is denoted as /a, b and it is defined as: /a, b = cos“^ ( ||a||j|b|| )’ where cos“^(-) stands 
for the inverse eosine funetion. For a random veetor a?, cc ~ CAf{i-i, S) denotes that x follows 
a Cireularly Symmetrie Complex Gaussian (CSCG) distribution with mean /x and eovarianee 
matrix S. Finally, Ei(-) represents the exponential integral function, which is defined in ll23] 
eq. (5.1.1)], while 7 ~ 0.5772 stands for the Euler-Mascheroni constant, as it is defined in ll23l 
eq. (4.1.32)]. 


II. System and Channel Model 

The spectrum sharing system, which is illustrated in Fig. [H is composed of a MISO incumbent 
system, comprising of a TX, TX 1, equipped with Mi antennas, along with its assigned single¬ 
antenna terminal, RX 1. Focusing on downlink communication, the incumbent system is willing 
to share its resources with a MISO licensee system. The latter system consists of a multiple 
antenna TX, TX 2, equipped with M 2 antennas, as well as of a licensee terminal, RX 2, assigned 
to TX 2. 

Considering the involved channels, spatially correlated Rayleigh fading is assumed for both 
direct and interfering channel links. As a consequence, for the channel between TX j and RX i, 
we have: hij ~ CA/’( 0 mi, R-iy)- 

The signal received at RX i, i ^ {1,2}, can be expressed as 

Vi = hfiWiSi + + Hi, ( 1 ) 

where, Wi denotes the transmit beamforming vector at TX i and it is assumed that Wi = 
y/PiUi, with Pi < and ll-UiH = 1 , where P^^^ is a maximum instantaneous power level 

at TX i. Also, Gaussian noise is considered at RX i, i.e., n* ~ CAA(0, Nq) and we assume 
that the information symbols for transmission are taken from a standard complex Gaussian 
codebook, i.e., s, ~ CAA(0,1), i e {1, 2 }. By analyzing ([I]), the instantaneous information rate 
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Fig. 1. The examined LSA system (post-licensing phase). 


of RX i, i e {1, 2} is given by [l2^ 


R^ 


loga 



]v„ + nih>5iV' 


( 2 ) 


In the seetion that follows, the problem of joint downlink preeoding with eombined, loeal CSI 
at the TX (CSIT), is formulated. 


III. Problem Formulation 

Foeusing on the deseribed system model, a realistie CSIT assumption that ean be made is 
that TX i, i G {1,2}, has both instantaneous and statistieal (eovarianee) knowledge of its direet 
links (i.e., TX 1 has instantaneous knowledge of direet link hi i and TX 2 has instantaneous 
knowledge of direet link ha, 2 ), whereas, the interferenee eross-links are merely statistieally known 
via knowledge of their eovarianee matriees. The seeond order statisties of the involved ehannels 
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constitute slow-varying information that can be realistically collected by each TX through low 
capacity/high delay links. 

Capitalizing on the available CSIT at TX i, i E {1,2}, the optimization problem of maxi¬ 
mizing the average rate of the licensee system, subject to an average rate constraint for RX 1 
can be formulated as di functional optimization problem, with functional dependencies related to 
the available CSI. Hence, the resulting optimization problem can be described as follows 

= argmaxE[i ?2 in 2 (^ 2 , 2 ))] 


subject to E [i?i m 2 (/i 2 , 2 ))] > Ti > 0, (PI) 

0 < 0 < ||U72(/12,2)||' < ^2““, 

where ri stands for the QoS demand of RX 1, in terms of average rate. 

The dependencies of the precoders to be optimized, on the corresponding instantaneous links, 
can be formulated by the following mappings 


Wi 


^ C 


Mi 


( 3 ) 


where i G {1,2}. Hence, the formulated optimization problem is of distributed nature. In what 
follows, for ease of exposition, we will omit to mention explicitly the dependencies of the 
precoders. 

One has to ensure that problem (IpTI) is feasible, in other words, QoS threshold ri has to be 
achievable when TX 1 transmits with full power MF precoding, in the absence of any interference 
coming from the licensee. This means satisfying the following 


E 


log2 1 + 




> n- 


( 4 ) 


The expectation over the interfering channels makes the optimization difficult to handle. However, 
exploiting the convexity of function log 2 (l -f 1), it becomes possible to apply Jensen’s inequality 
over the interfering channels. This significantly simplifies the optimization problem, while 





preserving its important features. The average rate expression for RX i, thus, beeomes 


E [Ri] — Eft,, 


log 2 1 + 


ml 




No + nmu. 




> E^.. 


loga 1 + 




No + Eh, j 


pim 


Ul 


(5) 


= E 




log 2 1 + 


No + PiV^Kau-i 


^ E 


Ri{wi,wi) 


Remark 1. It should be noted that this approaeh is only possible thanks to the faet that the 
preeoders Wi and W 2 are independent of the instantaneous eross-ehannels (as only the direet 
links are instantaneously known). □ 


With the aim of deriving a praetieal solution, slow power eontrol depending on the long term 
statistieal ehannel information, is assumed. Hence, instead of (instantaneous) power levels Pi 
and P 2 , we can use slow power allocation levels Pi and P 2 , where 0 < 

Altogether, in the remainder of the paper, we will work on the following optimization problem: 


(Pi*, ul, P 2 , U 2 ) = arg max E 


R 2 {Pi,Ui, P2, U2) 
subject to E Pi(Pi, lii, P 2 , 142 ) 

0 < Pi < 0 < P 2 < P^^^, 

||iii|| = 1, IIU 2 II = 1- 


> P, 


(P2) 


IV. Preliminary Results 

The following two propositions provide some characteristics of the optimal solution of prob¬ 
lem (lP2l) . 


Proposition 1. The ergodic rate constraint of RX 1 is satisfied with equality by any optimal 
solution of (lP2l) . i.e., 


E 


Ri{pfufp;,u*2) 


= p- 


( 6 ) 


Proof: The objective E P 2 (^i, P 2 , U 2 ) is monotonically decreasing with respect to Pi, 


while, on the other hand, the constraint E 


Pi (Pi, Ul, P 2 , U 2 ) is monotonically increasing and 

























continuous in Pi. As a result, one ean inerease the objeetive by redueing power level Pi up 
to the point, where the average rate eonstraint of RX 1 will be satisfied with equality. This is 
always feasible beeause ri > 0 implies that Pj* > 0. ■ 

The seeond proposition yields some insight with respeet to the optimal power alloeation 
seheme. 


Proposition 2. An optimal solution of problem (IP2I) satisfies that either TX 1 or TX 2 transmits 
with full power, i.e., when = p™'^^ or P 2 = 


Proof: Considering an optimal solution, one ean write P* = aiP, for some «! > 0 and 
P 2 = a^P, for some > 0, where P > 0. Then, taking every term of the objeetive and dividing 
the numerator and the denominator of its Signal to Interferenee plus Noise Ratio (SINR) by P, 
one obtains 


E 


R2{wI,W2) 


= E 


log2 1 + 




^ + at(n*fR2,int 


(7) 


whieh is a monotonieally inereasing funetion of P. Similarly, the aehievable average rate at 
RX 1 beeomes 


E 


Rl{wl, W*2 


= E 


loga 1 + 


a. 


\h^ 


u; 


f + «2l 


U 


H 


Ri,2'U2 


( 8 ) 


whieh is a monotonieally inereasing funetion of P, as well. 

If none of the two TXs transmits with full power, it means that it is possible to transmit 
with P' > P. Thus, the transmission using {a\P', uf Q^ 2 -P^ '“ 2 ) feasible and leads to a larger 
objeetive, whieh eontradiets the optimality of {alP, uf a^P, uf). ■ 


V. Statistically Coordinated Precoding 

We now present our main eontribution whieh is a new transmission seheme eonstituting a 
possible solution for optimization problem (iPTl) . Indeed, it is important to note that, although 
possibly suboptimal, our approaeh is able to guarantee the ineumbent rate eonstraint and is, 
therefore, a solution to the initial optimization problem. 


A. General Approach 

Sinee the derivation of elosed-form expressions for the optimal preeoders is hardly traetable 
due to the funetional nature of optimization problem (lP2l) (whieh requires optimizing over an 
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infinite dimensional space), we discretize the functional space and restrict the space of possible 
precoding solutions to a set of transmission strategies, S = {iSi,... ,5^}, where L denotes the 
number of different joint transmission schemes applied by TX 1 and TX 2. 

More specifically, each element. Si G S, refers to an ordered set P 2 ,Sii'^ 2 ,Si)^ ^ ^ 

Such a restriction to a finite set of joint transmission schemes allows for every transmit 
strategy. Si, / = 1,..., L, to be evaluated in terms of feasibility and in terms of performance. 
Furthermore, it provides a low complexity method for coordinating the TXs. 

1) Beamforming Design: The first step consists in designing the strategy set, i.e., the beam¬ 
forming strategies. Although any beamforming scheme could be chosen in theory, a good 
heuristic choice is key to the tractability and the efficiency of the approach. In this work, we 
restrict our analysis to the MF and the sZF strategies, as they represent the extreme approaches 
between which it will be necessary to strike a trade-off. 

MF precoding corresponds to the egoistic beamforming scheme, where TX i transmits using 


A 

Ui^MF = TTT-1 

r\ 


( 9 ) 


This beamformer maximizes the strength of the direct link without any consideration of the 
interference. 

In constrast, sZF corresponds to an altruistic beamforming scheme, where TX i transmits 
using 

sZF = ars max .^n. (10) 


The sZF beamforming scheme consists in exploiting the statistical information of the cross-links 
to reduce the created interference, while also taking into consideration the statistical information 
of the direct links. This strategy has the advantage of using only statistical information available 
at both TXs and hence enforces perfect coordination between the TXs, which will prove critical 
to an efficient joint transmission scheme. 

2) Power Control Policy: Power control is a key ingredient to ensure that the average rate 
constraint for the incumbent RX is not violated. Furthermore, it is shown in Section |W] that the 
incumbent QoS constraint is always fulfilled with equality and that one of the two TXs emits 
with full power, while the other reduces its power to respect the incumbent constraint. Therefore, 
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we denote by Vi the joint power policy where TX 1 emits with full power and by V 2 the joint 
power policy where TX 2 transmits with full power. 

3) Choice of the Transmission Policy: Considering the potential applicability of the two power 
control policies for each of the joint beamforming solutions, such a formulation leads to a 
joint transmission strategy set, which consists of 8 possible transmission schemes. However, the 
incumbent constraint is only fulfilled for some of the strategies and has to be verified otherwise. 
It is, hence, necessary to compute for each of these 8 transmission schemes the power emitted 
by one of the TXs and then evaluate the ergodic rate of both RXs. Once this is done, the best 
solution, in terms of average throughput for the licensee RX, is directly obtained. 

Remark 2. It is critical to understand that both TXs will always agree on which strategy to use 
as the TXs are statistically coordinated: only statistical information is necessary to evaluate the 
ergodic rates and choose the best strategy. □ 

B. Computation of the Ergodic Rates for each Strategy 

The ergodic rates for each of the 8 strategies need to be evaluated. However, the expressions 
are practically the same in the sense that the 8 possible strategies come from the combination 
of only a few parameters. We will hence only present in full detail two strategies: MF-MF-Pi 
and sZF-sZF-P 2 - The expressions for the other strategies can be trivially deduced. 

Remark 3. The feasibility of a given strategy has to be verified. However, the feasibility of 
the optimization problem is preserved as the feasibility is guaranteed for strategy MF-MF-Pi. 
Indeed, it contains the case where TX 1 transmits using MF and full power, while TX 2 does 
not transmit at all. □ 

1) Strategy MF-MF-Vi: The TXs transmit using the beamforming vectors Wi^mf and rt 2 ,MF- 
Furthermore, TX 1 transmits using Pi = It, thus, remains to determine how TX 2 controls 
its power to ensure that the incumbent ergodic rate constraint is fulfilled, i.e., that 


E[-Ri] > Ti- 


(11) 
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This can then be rewritten as 





( 12 ) 


where (a) holds by applying Jensen’s inequality to convex funetion log 2 (l + and the expee 


tation in the denominator ean then be eomputed using Lemma [3] in the Appendix with A = R 2_2 

and B = R| 2 R-i, 2 R-| 2 - 


Finally, a elosed form expression for the ergodic rate is obtained with Lemma [B Hence, the 


value of P 2 can be dedueed by biseetion, in order for the lower bound derived in (fT^ to be 
equal to ti. 

It remains to evaluate the eorresponding aehievable average rate of RX 2. Following a similar 
approaeh as the one for the ergodie rate of the incumbent, we ean obtain the following lower 
bound: 


E[R2] 



(13) 


Once more, the expectation in the denominator is obtained using Lemma [3l while a closed form 
expression for the ergodie rate is obtained with Lemma [IJ 

2) Strategy SZF-SZF-V 2 ' In this strategy, the TXs transmit using the beamformers Ui ^zf 
and 'U 2 ,szf, while TX 2 transmits using P 2 = It remains then to determine Pi- In that 
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setting, the rate of RX 1 can be lower bounded as 

PlIh.^l'Ui^sZFp 


E 




= E 


'^■1,1 


logs 1 + 


^0 + 


sZF 


> Ti. 


(14) 


This rate can be directly computed in closed form using Lemma l2l Finally, the power Pi, such 
that the ergodic rate constraint for RX 1 is met by the derived lower bound, can be obtained by 
bisection. 

It now remains to evaluate the corresponding ergodic rate of the licensee RX. This is given 
by the following expression 

Pr^\hl,U2,sZF\^ 


E 


R. 


= E 


h2,2 


logs 1 + 


No + Pi'u^,zpR2,iWi ,sZF 

The latter expression can be computed in closed form by applying Lemma |2l 


(15) 


VI. Reference Precoding Schemes 

In this section, we present two schemes which will be used to evaluate the efficiency of our 
statistically coordinated precoding approach. 

The first one, denoted as “interference temperature-based” precoding, is an adaptation of the 
approaches in the literature to allow for a fair comparison. Intuitively, it corresponds to the 
conventional underlay CR paradigm, where solely the secondary TX adapts its strategy in order 
for the interference received by the primary RX to be below a given threshold [jS]. 

The second one constitutes a coordination benchmark and it is a priori not reachable, but 
allows to bound the sub-optimality of the proposed approach. 


A. Interference Temperature-Based Precoding 

The interference temperature approach, extensively used in the CR literature, consists in 
forcing the secondary TX to create less interference to the primary user than a given interference 
threshold, which is here denoted by X. 

Considering that the secondary TX aims at minimizing the interference created and transmits 
using W 2 ,szf, the power emitted by the secondary TX is then given by 

ii D ’ ^2 

'“2,sZfR'1,2'U'2,sZF 


P 2 = min 


(16) 
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In order to conduct a fair comparison with the designed statistically coordinated precoding 
scheme, we need to determine the interference temperature, X, such that the ergodic rate con¬ 
straint of RX 1 is met with equality, i.e.. 


E 


logs 



pr^'Whi,! 

No+I 



= Ti. 


(17) 


The expectation appearing in (fTTI) can be computed by applying Lemma [T] The interference 
temperature threshold, X, can be then easily found by bisection. 


B. Coordination Benchmark 


When designing the beamformers, we can observe a clear trade-off between maximizing the 
desired signal (using MF) and minimizing the interference created. Hence, if we assume that 
one can achieve both goals at the same time, the following optimization problem is obtained for 
the power control, and leads to an, a priori, infeasible performance upperbound. 


max E logo 
Pl,P2 L 

subject to E 


1 + 


hwh 


2,2 I 


^0 + -Pi Amin (R-: 


2 , 1 , 


logs 1 + 


Plll^l,l 


No -f PsAmin (1^1,2; 


> X, 


(P3) 


0 < Pi < pr"^, 0 < Ps < Ps”^. 

The ergodic rate expressions appearing in (lP3l) can be computed in closed form by applying 
Lemma [TJ The optimal slow power control values are obtained by exploiting Proposition [2l 
Indeed, one of the two TXs transmits with full power, while the other one controls its power by 
bisection. Comparing the performance and the feasibility of both solutions leads to the solution 
of optimization problem (lP3l) . 


VII. Numerical Evaluation 

With the aim of evaluating the performance of the proposed statistically coordinated precoding 
scheme, extensive Monte Carlo simulations have been performed and, more specifically, 20000 
channel realizations have been simulated. We assume the existence of Mi = M 2 = M = 
4 antennas at each TX. Furthermore, we consider unit noise variance (Nq = 1) and a QoS 
threshold ri = 1 bps/Hz. 
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We consider a classical exponential channel correlation model [|25ll . in which the covariance 
matrices Rjj are given by 



1 

P 

P^ • 

1 

_1 




p 

1 

P 


, e {1,2}, 

(18) 


!>■ 

1 

pM-2 

pM-3 _ 

1 




where, ^ represents the pathloss and is chosen here equal to 1 when i = j and to 0.3 otherwise. 
In the investigated scenario the antenna correlation factor, p, is equal to 0.5. 



Fig. 2. Ergodic rate of RX 1 vs. transmit SNR, when incumbent QoS threshold n = Ibps/Hz. 

In Fig. 1^ and Fig. [3l the average rate of RX 1 and the average rate of RX 2 are depicted as 
a function of the system’s transmit Signal-to-Noise Ratio (SNR). The three curves represent the 
throughput performance achieved by the proposed statistically coordinated precoding scheme, the 
interference temperature-based precoding scheme, as well as the described coordination bench¬ 
mark. Focusing on RX 2, the coordination benchmark outperforms both the proposed precoding 
scheme as well as the interference temperature-based scheme, as expected. By observing Fig. [2l 
it should be noted that, in contrast with the coordination benchmark, the proposed precoding 
scheme fails to satisfy the incumbent average rate constraint with equality. This occurs because 
we resort to tackling optimization problem (lP2l) . which involves a lower bound of the average 
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Fig. 3. Ergodic rate of RX 2 vs. transmit SNR, when incumbent QoS threshold ri = Ibps/Hz. 


rate of RX 1. Nevertheless, the proposed algorithm suoeessfully manages to eontrol the average 
rate of RX 1 and this eapability ean be translated to a signifieant throughput gain for the lieensee, 
in eomparison to the one aehieved by the interferenee temperature-based preeoding seheme. 



Fig. 4. Ergodic rate of RX 1 vs. threshold n, SNR=10dB. 


The aehievable average rates of RX 1 and RX 2, by applying the proposed preeoding algorithm, 
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Fig. 5. Ergodic rate of RX 2 vs. threshold n, SNR=10dB. 

along with the ones aehieved by the two referenee preeoding sehemes, are depleted in Fig. [Hand 
Fig. m respeetively, as a funetion of QoS threshold, ti, when the transmit SNR of the system 
is equal to 10 dB. The average rate eonstraint for RX 1 is fulfilled by all three sehemes for the 
whole examined range of ri. Also, the proposed preeoding seheme outperforms the interferenee 
temperature-based one, espeeially when the average rate eonstraint of the ineumbent is loose, 
whieh oeeurs due to the faet that under this regime there is more to gain for the lieensee by 
means of an effieient eoordination. 


VIII. Conclusions 

In this work, we propose a novel, joint preeoding seheme with referenee to a shared speetrum 
aeeess system where the two TXs eoordinate on the basis of statistieal knowledge of the global 
multi-user ehannel. Our approaeh eonsists in formulating a Team Deeision problem, the solution 
of whieh is approaehed by redueing the transmission strategy spaee to a finite number of 
strategies. This method is key to enforeing between the TXs and obtaining a praetieal solution to 
the intrieate Team Deeision problem. Sueh an approaeh allows to improve over the eonventional 
underlay CR approaeh, by enforeing more eoordination between the two eoexisting TXs at 
the priee of low CSI and eomplexity requirements, as the eoordination ean be realized offline. 
Approaehing the global optimum is both a diffieult and ehallenging problem that ean be taekled 
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in the future. The proposed seheme has also a strong potential in other more eomplex seenarios 
with multiple ineumbent and/or lieensee networks. 


Appendix 

Lemma 1. I^S\ eq. (37)] Let 7 G M"*" and h ~ CA/'(0n,Rh), where covariance matrix R^, has 
n distinct eigenvalues It then holds 


Eh [log 2 (l +711^11^)] = 


E 




( 19 ) 


in(2)7n”.Di^Gn:.wdi-A)' 

Lemma 2. KTTl eq. (75)-(76)] Let 7 G M"*" and w G be deterministic, and h ~ CJ\f{0n, Rh,)- 


It then holds 


E, 




loga I + tI^ 


ln( 2 ) 




7 Al(Re#)/ ’ 


( 20 ) 


where Ai(Re^) is the unique non-zero (positive) eigenvalue of matrix Rgj = R^mm^R^. 


Lemma 3. Let us consider two positive semi-definite matrices A and B in with eigenvalues 
denoted as Ai(A),..., A„(A) and Ai(B),..., A„(B), respectively, where it is assumed that A is 
of full rank and has no multiple eigenvalues. We also assume that matrix A can be decomposed 
as A = UaAaUH, where Ua A a unitary matrix and Aa = diag (Ai(A),..., An(A)), where 
it holds that 0 < Ai(A) < ... < A„(A). Let x G be a standard complex Gaussian random 
vector, such that x ~ CAA(0„,I„). It then holds 

Ac^Bccl I Ai(A)'"-2 _ ;L) (ln(Ai(A)) - 7 ) + 1 ) 


E 


x^Ax 


E 


i=l 


(Ai(A) — Aj(A)) 


A,(A)--i (ln(A,(A)) - 7) (A*(A) - A,(A)) 


E 

k=l,k^i 


(^*(A) — Aj(A))^ 

A,(A)-i (ln(A,(A)) - 7) ih{A) - A,(A)) 


( 21 ) 


HjjLk (Afc(A) - Aj(A))j 


where B = U^BUa- 

Proof: We prove this result in two steps. Firstly, we show that considering two matrices A 
and B with different eigenbases, we can come back to the case of matrices having the same 
eigenbasis. We then prove the lemma for this case. 
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Let us assume that A and B have different eigenbases. We consider their eigendecompositions 
A = UaAaU^ and B = UeABUg, where the diagonal entries of A a are sorted in an 
increasing order and the diagonal entries o/A b are non-decreasingly sorted. Introducing matrix 
B = U^BUa, the expectation in question becomes 


E 


x^^x 

— F 

Ac^U^BUAai 

— F 

x^^x 

— F 

\Y.t.Y.U[B\,xtx,-\ 

x^Ax 

11 

x^A^x 


x^Aa.x 


ELi Afc(A)|xfc|2 


( 22 ) 


Exploiting the fact that x ~ CA/'(0„,I„), if we write each x*, i = 1,... ,n in polar repre¬ 
sentation, i.e., Xi = then we have that all phases fi, i = 1,... ,n and amplitudes 

\xi\, i = 1,... ,n are mutually independent and the phases are uniformly distributed. As a 
result, the expectation takes the following form 


E 


x^^x 

— F 


— F 

\ 1 

x^Ax 


ELi Afe(A)|xfc|2 

II J 

ELi Afc(A)|xfc|2 


(23) 


or, equivalently 


E 


x^^x 

= E 

x^JJa diag 

[B] 

11 ,..., 

1 

ffiC 

x^Ax 


x^Ax 



(24) 


Hence, the case of equal eigenbases is recovered. 

Consequently, we proceed by considering without loss of generality that matrices A and B 
have the same eigenbases. However, it should be noted that elements i = 1,... ,n are 

not sorted in any particular order. 

Focusing, now, on the derivation of a closed form expression of the expectation, we have that 


x^Ax 


E 


= E 


Er=i 

E-=iA,(A)|x,f 

d 


i=l 

n 

E 1^1.,. 


SAi(A) 

a 


E [-bIm® 

i=l 

n 


Xi 


E-=iA,(A)|x, 


2=1 


aA,(A) 


E 


0=1 

n 


(25) 


i=i 


Xi 


Let us define random variable (RV) X = J2j=i Aj(A)|xjp. Using ^2^ eq. (8)], it is shown by 
induction that the Probability Density Function (PDF) of X is the following 

P^W = En - 


k=l j^k 


Afc(A)-Aj(A) 


(26) 
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As a result, the expectation of RV ln(X) is given by the expression that follows 

Afc(A)"-2 


poo 

E [ln(X)] = ^ dx'^^ 

k=i j^k 


A.(A)-A,(A)- 


Exploiting /[^ eq. (4.331.1)], expression (1271) becomes 

Afc(Af-i(ln(Afc(A))-7) 


E |ln(.Y)] = 5^ ^ 

k=l 

Taking, now, the partial derivative of 


d d 

E|lnmi = 


(A.(A) - Aj(A)) 

, with respect to A*(A), we obtain 

A,(Ar-i (ln(A,(A))-7) 


aA,(A) 


5A,(A) 1 n,v.(A.(A)-A,(A)) 


d 


d\{A) 


E 


Afc(Ar-i(ln(Afc(A))-7) 


For the first term of 


d A,(Af-i(ln(A,(A))- 7 ) 


k=i,k^i njVfc(Afc(A) Aj(A)) 

I, the following expression is obtained 

Ai(A )^-2 _ X) (ln(Ai(A)) - 7 ) + 1) 


5A,(A) 1 n,v.(A.(A)-A,(A)) 


(27) 


(28) 


(29) 


Uj^i (Ai(A) — Aj(A)) 

A,(A)--^ (ln(A,(A)) - 7) “ A, (A)) 


n,^,(A.(A)-A,(A)) 


(30) 


It now remains to find the second term of (l29l) in closed form. We, thus, obtain the following 

d \ ^ A,(A)"-i(ln(Afc(A))-7)l A d [ Afe(A)’^-i (ln(A,(A)) - 7 ) | 
n,,JA»(A) - A,(A)) / ^2^^7A.(A)\ n,^.(A.(A)-Aj(A)) / 

(31) 

Given that i f k, the partial derivative appearing in the right hand side of (l3T]) . is given by the 
following expression 


d J 

fAfe(A)"-i(ln(Afe(A))-7)l 

I Afe(A)"-i(ln(Afc(A))-7)^j 

(Afc(A) — Aj(A))| 

^X^{A.) 1 

[ rij/fe (Afc(A) — Aj(A)) J 

f 1 

(Afc(A) — Aj(A))^ 

^2 


(32) 


qa^a) (Afc(A) - Aj(A))| = - (Afc(A) - Aj(A)). This concludes the proof 
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